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A general, exact formula is derived for the expectation value of the electromagnetic energy density 
of an inhomogeneous absorbing and dispersive dielectric medium in thermal equilibrium, assuming 
that the medium is well approximated as a continuum. From this formula we obtain the formal 
expression for the Casimir force density. Unlike most previous approaches to Casimir effects in 
which absorption is either ignored or admitted implicitly through the required analytic properties 
of the permittivity, we include dissipation explicitly via the coupling of each dipole oscillator of the 
medium to a reservoir of harmonic oscillators. We obtain the energy density and the Casimir force 
density as a consequence of the van der Waals interactions of the oscillators and also from Poynting's 
theorem. 

PACS numbers: 42.50.-p, 03.70.+k, 



I. INTRODUCTION 

Analyses of quantum fluctuations of the electromag- 
netic field in dielectric media typically take a macroscopic 
approach in which the field is quantized under the as- 
sumption that the medium is characterized by an elec- 
tric permittivity [l[. Much of this work is restricted to 
fields at frequencies for which absorption is negligible, 
and some of it assumes furthermore that dispersion is 
negligible. An important exception is the work of Hut- 
tncr and Barnett [2|, who allow for both absorption and 
dispersion; as in virtually all work in this area, they 
treat the (unexcitcd) atoms of a homogeneous dielectric 
medium |3| as harmonic oscillators [J] , and dissipation is 
accounted for by coupling each of these oscillators to a 
"bath" of reservoir oscillators. The atoms are also cou- 
pled to the electromagnetic field, and the entire system 
of atoms, reservoirs, and field oscillators is then diag- 
onalized along the lines of Fano's method [J to yield 
expressions for quantized electric and magnetic fields. 

In a less elegant approach one writes Heisenberg equa- 
tions of motion for the atom, reservoir, and field vari- 
ables. The effect on each atom of its coupling to its reser- 
voir is to introduce a damping force and a Langcvin force, 
the latter ensuring the preservation of canonical commu- 
tation relations. The Langevin force on each atom results 
in a fluctuating "noise polarization" of the type intro- 
duced in Rytov's theory of fluctuational electrodynamics 
Q , and used by Lifshitz Q in his theory of the van der 
Waals-Casimir force per unit area between two dielec- 
tric half-spaces. The operator Maxwell equations with 
this noise polarization yield expressions for the quantized 
transverse fields having the same form as those obtained 
by Huttner and Barnett 0,11]. 

One purpose of the present paper is to extend earlier 
work i^, lEj , hereafter referred to as I, to the case of inho- 
mogeneous dielectric media. This leads to a formula for 
the total energy density of a dispersive, dissipative, and 



inhomogeneous dielectric medium in terms of the dyadic 
Green function G. We derive the force density at finite 
temperature T associated with spatial variations of the 
complex electric permittivity e(r, a;) of the medium: 
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Gii{Y,Y,Lo) is the sum of the three diagonal components 
of the Fourier transform of the Green dyadic in the sum- 
mation convention for repeated indices used throughout 
this paper. From ([T]) one obtains, for instance, the Lif- 
shitz formula [6| for the force per unit area between two 
semi-infinite dielectric media separated by a distance d, 
as was shown in the earl y w ork of Dzyaloshinskii et al. 
\]S\ and Schwinger et al. [11| , and from that formula one 
obtains in various limits, as is well known, the Casimir 
force between perfectly conducting plates, the van der 
Waals interaction between two atoms, the Casimir-Polder 
interaction between an atom and a conducting plate, etc. 
The main motivation for the present work is not just to 
rcdcrive the general force density ([1]) by different meth- 
ods, but rather to obtain it by allowing from the outset 
for absorption. Derivations of the (T = 0) Lifshitz for- 
mula based on changes in electromagnetic energy J12l |. 
following Casimir's original calculation [IJ], are based on 
the zero-point electromagnetic energy ^ ■ hiUj/2, the sum 
being over all possible modes of the field. The frequencies 
LUj in such calculations are first determined for the case 
of nondissipative media (real permittivities), and dissi- 
pation is later accounted for by making in effect an ana- 
lytic continuation based on the requirement from causal- 
ity that the permittivity is analytic in the upper half of 
the complex frequency plane. The question of how to ap- 
ply this approach by admitting dissipation from the start 
is an old one in the theory of Casimir forces. Agarwal 
J14l | , for instance, remarks that "if the damping of the di- 
electric function is included, then the [method based on 



zero-point energy] seems to fail. In presence of damping 
the normal-mode frequencies are complex, and it is not 
clear what one should sum over to obtain the interaction 
energy." In Reference [l^l, similarly, it is noted that "it 
is not obvious how to extend [the zero-point-energy ap- 
proach] to the case of absorbing media, where the [per- 
mittivities] arc complex." Ginzburg jl6| observes that, 
"Apart from a general maxim that 'victors should not 
be judged' one can justify the derivation of [the Lifshitz 
formula] on the basis of [the assumption that a mode 
of frequency lo has zero-point energy hiu /2] for absorbing 
media as follows. Firstly, the permittivities [appearing in 
the Lifshitz formula] are functions. Secondly, the func- 
tion e(w) is always real on the imaginary axis. The result 
obtained for transparent media ... must therefore clearly 
be the same as the appreciably more general one which 
is applicable to absorbing media. However, this would 
hardly confirm such a conclusion unless it had been ob- 
tained earlier without additional assumptions. Both for 
this reason and also bearing in mind other related or 
similar problems one must somehow consistently gener- 
alize the expansion in cigenoscillations with frequencies 
LOa ... to absorbing media." Ginzburg proceeds to intro- 
duce the concept of (orthogonal) "auxiliary" field modes 
associated with the immersion of the entire system "in 
some auxiliary resonator with perfectly conducting walls 
... The frequency w is considered to be a parameter, while 
the eigenfrequencies of the resonator Wq are determined 
from the homogeneous field equations" corresponding to 
Maxwell's equations without a fluctuating polarization 
density. 

There are arguably more direct ways of explaining the 
success of calculations that "consider directly only trans- 
parent media" and attribute Casimir forces to changes in 
zero-point energy due to the presence of material media. 
A clue in this direction is provided by the case of a homo- 
geneous dielectric medium: the zero-point energy has the 
same form regardless of whether one allows "directly" for 
absorption [8[. 

In the following section we derive the Gasimir free en- 
ergy and force density for a dispersive and absorbing, 
linear, inhomogeneous dielectric medium based on an ex- 
tension of the classical formula 



the seminal work of Dzyaloshinskii et al. [T^j using dia- 
grammatic methods and subsequently by other authors 
by various other techniques. In Section [ITT] we rederive, 
from Poynting's theorem, the expression for the energy 
density obtained in Section |TT1 and show that it reduces 
to the result obtained in I in the special case of a homoge- 
neous medium in which the complex permittivity has no 
spatial dependence. Section lTVl prcsents a brief derivation 
of the Casimir force density for a dispersive and dissipa- 
tive dielectric based on the Maxwell stress tensor. Our 
conclusions are briefly summarized and discussed further 
in Section |Vl As in I we focus mainly on the case of zero 
temperature. The extension to finite equilibrium temper- 
atures is straightforward and so we simply summarize it 
in the Appendix, where we also briefly outline the deriva- 
tion of the Lifshitz formula based on the force density ([T]) 
or its finite-temperature generalization. 



II. ENERGY AND FORCE DENSITY FROM 
COUPLED DIPOLES 

A. Interaction Energy in Terms of Free-Space 
Green Dyadic 

We start by considering a collection of M electric 
dipoles in free space. The nth dipole ( "atom" ) has associ- 
ated with it an electric dipole moment operator p„ = ex„ 
and a resonance frequency wq , and it is coupled to a reser- 
voir that results in a damping rate 7 and a Langcvin force 
Fi,i. It is also coupled to the total electric field E(r„,t) 
at its position r„. The Heisenberg equation of motion 
for the electron coordinate operator for the nth atom is 
derived in I [Eq. (51)] [l^ 



1 - e - 

Xni +7*m +Wo*m = —FLni{t) -\ £'i(r„,t). (3) 

m m 



.2 - 



The first subscript in) on x identifies the nth atom, while 
the second subscript (i) denotes the ith Cartesian com- 
ponent of x„. As in I we use a circumficx to denote 
quantum-mechanical operators, and we define Fourier- 
transformed operators p„i(a;) = ex„i(aj), etc. by writing 



W, = -ia(w)E2(r) 



(2) 
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for the change in energy when an electric field E of fre- 
quency u) induces an electric dipole moment a(w)E in a 
particle having a real polarizability a{uj). For an atom 
in a state characterized by the polarizability a{u)) this 
is just the quadratic Stark shift of its energy level. The 
generalization needed to obtain the Gasimir free energy 
only requires allowance for a complex polarizability and 
the effect on each dipole of the fields from all the other 
dipoles. We explicitly include dissipation and Langevin 
forces resulting from the coupling of the dipole oscilla- 
tors to their reservoirs and show that the force density 
obtained in this way is equivalent to that obtained in 



FLm{t)^ dLj[FLm{L0)e-''^' + FLm{~Lu)e"'% (5) 
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with p„i(— w) — pI^^{uj), etc., which follows from the re- 
quirement that the operators p, Fl , and E be Hermitian. 
From ©, 

[wg -Up- - i^ij:\pm{uj) = — i^L„i(a;) H Ei{Vn,uj). (7) 



The electric field in (O is the sum of the source-free 
("vacuum") field Eoi{rn,i^) at r„ and the fields from all 
the dipolc sources: 

Ei{rn,uj) = £'oi(r„,w) + — ^ G'°,-(r„, r™, w)p„y(a;), 



(8) 
where G°(r, r',aj) is the Fourier transform of the (re- 
tarded) free-space dyadic Green function satisfying 

V X V X G°(r, r') - ^G°(r, r') = AnS^r - r'). (9) 
Equation ([7]) is therefore 
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or, in matrix form, 

p(a;) = A{uj) 



-Fl{cj) + -Eoiuj) 
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where pioj), Fl(uj), and Eo{uj) are 37\A-dimensional col- 
umn vectors {J\f dipoles, 3 Cartesian coordinates) and 
the 3M X 3A^ matrix 



A{u:) 
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We are interested in the change in energy involved in 
bringing the JV dipoles from a configuration where they 
are infinitely far apart and not interacting to one where 
they are separated by finite distances and interacting 
with themselves as well as with the source- free, fluctu- 
ating electric field. Consider first the simple example of 
a single dipole with real polarizability a{uj). According 
to the Hellmann-Fcynman theorem [18| the change in en- 
ergy when the dipole is brought from infinity (E = 0) to 
a point r where the electric field of frequency w is E is 



W = - 



^ ^(p.Et) = -^ ^«(c.,A)(E(r).Et(r)), 

. . (^^) 
where the coupling constant in the integrand, in this case 

the electric charge, is taken to be Ae. Since a(w) is pro- 
portional to e^, i.e., a{uj,X) = A^a(w), Eq. ^^ simply 
generalizes the classical formula ^. The relation [I3 
[see also Eq. ([57]) below] between the vacuum expecta- 
tion value of the free-field operator Eo(r, oj) • Ej(r,a;') 
and the Green dyadic Go 

(Eo(r,w)-ES(r,c^')) - -ImTrGo(r, r,c.)<5(w-w'), (15) 



where Tr is the 3x3 trace, suggests that 

W= / ^ / «ia;a(a;,A)ImTrGo(r,r,( 
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when the field is the source-free electric field including 
all frequencies. But of course the polarizability cannot 
be real at all frequencies, and the correct form of W is 
not (HSI but [H 

W = Im / — / da;Tra(a;,A)Go(r,r,w). (17) 
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Thus if we again use a(a;,A) = A^a(aj), we recover the 
expression 
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dwTra(w)Go(r,r,w) (18) 



obtained by other methods [20| . 

For our collection of dipoles we identify from (fTTj) the 
polarizability matrix (e^ lm)A(uj) and define the expec- 
tation value of the interaction energy as [21| 
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where we define the polarizability 



aoiuj) 



UJq ~ Lo"^ — ijuj 



(19) 
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(Spatial coordinates r„,rm are treated here along with 
Cartesian components i, j as matrix indices.) 

The integrand in (jl9|) can be expanded in powers of 
ao: 
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The first term in brackets is part of a single-particle self- 
energy, while the terms that arc non-diagonal in the space 
coordinates r„ , r^ correspond successively to two-body, 
three-body, etc. van der Waals interactions [10, H, [22, 
[23J . Thus, for instance, the second term in brackets, 
written out explicitly using 



1 1 -^ 

TrlgO^o^oGo] = 2 X! Q;o(^)Gojj(rn, fm, w) 

771, n—1 
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is found to be just the sum of pairwise van der Waals 
interaction energies of the TV atoms when terms with 



m ~ n are excluded (these correspond to self-energies). 
In the model in which the atoms form a continuum we 
replace the summation in (|22|) by 



o]^ (23) 
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where N{v) is the number of atoms per unit volume at 
r. Using the formula 



e{uj) = 1 + ^TiNaaiuj) = 1 + 
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for the permittivity, we similarly replace (|2ip by 

"e-1" 
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(25) 



(26) 



in the continuum approximation [24| . We have presumed 
in deriving this expression that the polarizability ao does 
not depend on r, so that the dependence of the permittiv- 
ity on r stems solely from the r-dependence of the number 
density N{r). However, it is straightforward to rederive 
([26|) with an r-dependent ao, so that different parts of 
the medium can have different resonance frequencies as 
well as different number densities. 



B. Energy Density in Terms of the Full Green 
Dyadic 

G° satisfies ([9]), while the Green dyadic G in the case 
of a medium with complex permittivity e(r,a;) satisfies 



V X V X G(r,r') 



-G(r,r') = 47r53(r-r'), (27) 



together with appropriate boundary conditions. There- 
fore 

Vx Vx (G-G°)-^(G-G°) = ^[e-1]G, (28) 

and the solution of this equation obtained using the 
Green function G° implies the Born-Dyson-type relation 



G^G° + G° 
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I.e. 



G„-(r,r,c.) - G°,.(r,r,o.) + ^y'dVG°,(r,r',c.) 
„ e(r',w)-l 
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From iMl) and IMl) it follows that 



u{r) = -—Im dc.Trlog[(G")-iG], (31) 
27r Jo 

which is the well-known "trace- log formula" [25[. The 
free-space Green dyadic G'^ is independent of the atoms 
or the properties of the medium formed by them. The 
term log[(G°)^^] in the above expression therefore sub- 
tracts from the total Green function G the "bulk" con- 
tribution G", and therefore u(r) contains only the "scat- 
tering" part of the Green function, that comprises the 
interaction between the atoms. 

We can express this energy in a different form by par- 
tial integration as in Reference |2al . Omitting the "bulk" 
contribution in (|5T|) . 



ulv) ~ — ImTr / dwo;— — logG 
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(32) 



where we use a prime to denote differentiation with re- 
spect to uj. Now from (|27| . 



(33) 



^2 ^2 

V x V x ,5G - -^eSG = -^SeG, 
c c 



and therefore, formally. 
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where of course there is an integration over space implied 
on the right-hand side. Similarly 
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V X V X G' - — eG' = ^eG + ^e'G, (35) 






and therefore 



Using this result in (|32p . we obtain 
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(37) 



This is our general expression for the energy density. We 
discuss it further and give an alternative derivation of it 
in Section [nil 

Two points concerning the derivation of p7|) arc worth 
noting here. First, we have seen that the energy density 
pG]) is associated with many-body van der Waals inter- 
actions, and therefore vanishes when the polarizability 
Q!o(a;) — > and the permittivity e{uj) — >■ 1 for all fre- 
quencies, as can also be seen from (|3ip . However, when 



e{ij) — ^ 1 the energy density 
(PS)) becomes 



37l) that was derived from 
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where we have used the fact that TrImG'°(r, r, uj) = 2uj/c. 
Thus we recover the correct vacuum electromagnetic en- 
ergy density in the hmit e{uj) — >■ 1. The reason for this 
is that in obtaining the trace-log formula (|32|) from ([3T|) 
we have effectively added the part of the energy den- 
sity needed to make u(r) have the correct (nonvanishing) 
vacuum- field energy density when e(cj) -^ 1. 

The second point is that there are other contributions 
to the total energy density that are not included in ((37)) 
and that do not contribute to the Casimir force density. 
Of course this is not surprising, as in deriving (|37|) we 
considered only the energy involved in the induction of 
dipole moments by the electric field and in the interac- 
tion of these dipoles. Thus our derivation does not ac- 
count for the constant energy (density) absorption rate 
of a medium with permittivity e(r,a;); this absorption 
rate -Rabs has exactly the form expected from classical 
electromagnetic theory: 



R 
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-Tr / dwwe/(r,w)(E(r,w)-Et(r,w)), 
-^-^Tr / dwiJ^eiir,Lo)lmG{r,r,Lo), (39) 



where in the second line we have used the generalization 
of (fT5|) given in (|87p below. This fiow of power into the 
atoms' reservoirs is of course exactly cancelled by the 
power lost by the field, so that there is no net change in 
total energy density. In the calculation using the Poynt- 
ing theorem in the following section, all contributions to 
the total energy density are included, and this cancella- 
tion is seen explicitly. 



In addition to the fluctuating dipoles induced by the 
electromagnetic field, whose interaction is responsible for 
van der Waals (Casimir) forces, there are fluctuating 
dipoles due to the Langevin noise forces acting on the in- 
dividual atoms. These "noise" dipole moments do not re- 
sult in forces among the atoms, since the Langevin forces 
acting on different atoms, unlike the electric fields induc- 
ing dipole moments in different atoms, have no spatial 
correlations. In the presence of the fluctuating electric 
fleld they do, however, contribute (K = 47rPnoise is the 



noise polarization deflned in Eq. 
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(40) 



to the total energy density. This is a single-particle 
self-energy resulting in effect from the interaction of the 
Langevin-force-induced dipole moment of each particle 
with the part of the electric field due to this same dipole 
moment. As such it has the effect in our model of deter- 
mining in part the equilibrium positions of the particles 
within the medium, but does not contribute to interpar- 
ticle interactions. 



C. Force Density 



The force density can be obtained from either expres- 
sion (|26)) or (|37)) for the energy density. Consider first 
(P51) . The fact that G" is independent of any properties 
of the medium means that a variation S£ of S due to a 
deformation of the medium depends only on the variation 



Se of e that accompanies the deformation jT3] ■ Therefore 
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It follows from pop that the quantity in curly brackets is 
just the trace of the Green dyadic G(r,r,w), and conse- 
quently the variation in the energy density is 



6u{r) 



S-K^c^ 



Im / da;w2(5e(r,w)TrG(r,r,a;). (42) 
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Following Dzyaloshinskii et al. [lO[ , we consider an in- 
finitesimal local displacement R(r) of the medium. This 
transport of the medium implies a variation in e at r such 
that e(r,cj)-|-5e(r, w) = e(r — R, w), or (5e(r, w) = — Ve-R. 
Therefore 



Su(r) = --— ^^Im / da;w2[-Ve-R]TrG(r,r,a;) 



d'Vf (r) • R, 



(43) 



where the Casimir force density f (r) obtained in this way 
is given by equation ([T]): 
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where we used that Tr G = Gu . This is equivalent to 
the result of Dzyaloshinskii et al. obtained by summing 
a sequence of diagrams, each successive one including a 
number n of closed loops corresponding to n-body van 
der Waals interaction energies. As these authors discuss, 
equation ([T|) does not in general give the total force den- 
sity, as one must account for the variation in e due to 
changes in density as well as to the displacements con- 
sidered in obtaining (j43p . Ignoring the former amounts 
to assuming displacements R such that V • R = 0. For 
the calculation of forces on bodies embedded in uni- 
form fluids in mechanical equilibrium, however, the addi- 
tional contribution to the total force resulting from den- 
sity variations must be balanced by pressure arising from 
non-electromagnetic forces, and may in effect be omitted 
jl(j |. It must also be noted that G(r,r,a;) diverges un- 
less spatial frequencies are cut off at some large value 
on the order of 1/ (interatomic spacing). However, again 
as discussed by Dzyaloshinskii et al.., the contributions 



from high spatial frequencies are the same at each point 
r for an inhomogeneous medium as for a homogeneous 
medium having the same value of e at r as the inhomo- 
geneous medium. The divergent Green dyadic appearing 
in ((!]) should therefore be replaced by [lO| 



G(r,r,w)-G(r,r,Lj), 



(45) 



where G is the Green function of a homogeneous medium 
with the same value of e at r as the inhomogeneous 
medium under consideration [27|. This subtraction of 
the singular G(r, r, w) is discussed further below. 

Consider now the derivation of the force density from 
the energy density as expressed by ([57]) . A variation Se 
in the permittivity results in a variation 



6u{r) = ^-^-^Im / dujLu^\[2Se + ujSe']Gi., 
+ [2e + 0Je']5G^^} 
in u(r). From ^ and ((36)) . 
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[2Se + ujSe']G + [2e + uj€']SG = [2Se + LuSe']G 

1 LO^ 

+ \2e + uje'] ^,5eGG = \2Se + ojSe']G + ujSeG' 

= 26eG + Lj^[5eG]. (47) 
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Partial integration of (|46|) then gives 



6u{r) — — 



Stt^c 



-Im 
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implying again the force density ([Ij when e varies with 
position. Essentially the same calculation was presented 
by Milton ct al. [23 under the assumption that e/ = 0. 
An "inhoiTLOgeneous medium" as defined here includes 
the important case of spatially separated homogeneous 
bodies, as in the case of the two parallel dielectric slabs 
treated by Lifshitz 0. Except for such simple geome- 
tries, the calculation of Casimir forces based on ((U — 
which amounts in effect to the calculation of the (classi- 
cal) dyadic Green function — must be performed numeri- 
cally. Our interest here, as discussed in the Introduction, 
is only in the general question regarding what effect dissi- 
pation has on the derivation of ([T]) by different methods, 
especially in derivations based on energy variations. 



D. Remarks Relating to Early Work 

In their diagrammatic approach Dzyaloshinskii et al. 
obtained the following expression for the force density (in 
addition to the aforementioned contribution that arises 
from displacements R with V • R 7^ 0): 



f(r) 



ksT 



^w2Gi,(r,r,iw„)Ve(r,ia;„). (49) 



Here the (Matsubara) frequencies aj„ = 2TrnkBT/h, and 
the prime on the zero below the summation sign means 
that a factor 1/2 multipHes the n = term. For T — !> 
the sum can be replaced by an integral, and P^ becomes 



f(r) 



87r2c2 



dee'Ve(r,ze)G,,(r,r,je). (50) 



G, like e, is analytic in the upper half of the complex fre- 
quency plane. Using this fact, it is easy to show that ([50)1 
is equivalent to ^ . The appearance of the permittivity 
at imaginary frequencies in the work of Dzyaloshinskii et 
al. can be traced to an analytic continuation of a Green 
function obtained with e evaluated at real frequencies 
to one in which e is evaluated at imaginary frequencies. 
Thus, consistent with the remarks in the Introduction 
about e being a function that is real- valued on the imag- 
inary axis, the validity of the derived force density for 
dissipative media rests on the formal analytic properties 
of e. The same is true of the analysis of Schwinger et 
al. [11|. Our more "brute- force" approach does not ex- 
plicitly invoke these properties; we calculate the Casimir 
free energy by exhibiting explicitly, albeit with a simple 
model, the dissipative character of the dielectric. 

Mahan [221 considered essentially equation ([3|) without 
damping and Langevin forces and without coupling to 
the fluctuating source-free field. The eigenfrequencies ftj 
of the coupled-oscillator system are obtained by writing 
x„ = x„(0) exp(— zfit) and then applying the argument 
theorem to obtain the difference in the zero-point ener- 
gies X]i(l/2)^^j between a given configuration and that 
when all the interatomic distances are infinite. Renne 
[23| noted inconsistencies in this approach in the case of 
retarded interactions, and proceeded in a similar fashion 
by solving for the eigenfrequencies of the coupled system 
of dipole oscillators and field oscillators. 

Renne includes radiative reaction, as does Agarwal 
[ij] in a general response-function formulation, and con- 
sequently the polarizabilities of the material oscillators 
arc complex; radiative reaction makes the only contribu- 
tion to the imaginary part of the polarizability in their 
work. However, this imaginary part is connected not 
with dissipation of field energy as such but with (energy- 
conserving) Rayleigh scattering of radiation, consistent 
with the optical theorem. To compare with the result of 
Renne, for instance, for the interaction energy of a sys- 
tem of ground-state harmonic dipole oscillators, we use 
the identity Tr log X = log detX IHl to write ^ as 

t /'OO 

£ = -— Im / da;logdet[l — ao('^)Go(w)] 
27r Jo 

= — / d^logdet[l-ao(*e)Go(*0], (51) 

which is equivalent to equation (18) of Renne [23[. In 
his derivation, however, the damping rate 7 in the defi- 
nition (|20p of ao(a;) is obtained from an approximation 
to the radiative reaction field in which the third deriva- 
tion with respect to time of a dipole moment is replaced 



by minus the square of an oscillation frequency times a 
first derivative with respect to time, and in this approx- 
imation the imaginary part of q:o(w) provides in effect 
for a coUisional- or ohmic-type dissipation rather than 
for elastic scattering of radiation. (Without this approx- 
imation the permittivity in Renne's approach would not 
be analytic in the upper half of the complex frequency 
plane, and it would preclude the analytic continuation 
of a Green function that was alluded to earlier.) In our 
approach 7 is obtained, without approximation, from the 
coupling of each oscillator to its reservoir, and e(a;) [Eq. 
(1211)] has no contribution from radiative reaction. Dis- 
sipation of field energy is due, as usual, simply to the 
transfer of energy from the field to the atoms and the 
subsequent transfer of this energy to the reservoirs rather 
than back to the field. 

Radiative reaction is of course naturally included in our 
approach; it is associated with the term G°(r„, r„,a;) in 
equation (jlOp and with a rcnormalizable self-energy. The 
complex permittivity p4p appearing in equation ()27|) for 
the Green function, however, has no contribution from 
radiative reaction. This is related to the implicit assump- 
tion in our continuum model that there are no density 
fluctuations and therefore that the extinction coefficient 
due to Rayleigh scattering vanishes [30|. In the case of 
two spatially separated homogeneous bodies, for exam- 
ple, the only "scattering" that occurs is in the form of 
reflection and refraction at boundaries. 



III. MACROSCOPIC QED APPROACH TO THE 

ENERGY DENSITY OF A DISSIPATIVE 

DIELECTRIC MEDIUM 



In I we derived expressions for quantized electric and 
magnetic flelds and for the energy density in a dispersive 
and absorbing, homogeneous dielectric medium. We now 
extend these considerations to inhomogeneous media. 



A. Poynting's Theorem 

As in I we base our approach essentially on the macro- 
scopic, Heisenberg-picture Maxwell equations and the 
Poynting's theorem that follows from them. In the 
conventional notation Poynting's theorem for the sym- 
metrized Poynting vector operator S = (c/87r)[E x H — 
H X E] takes the form 



(S)-nda 



hb^ 



/ H 



<9D 

dt 


+ 


dt 


dt 


+ 


dt 



E)d-^r 



ii)d^r 



-jjui..t)d\ 



(52) 



or, after a time-integration, 



where 



dt' f (S) -nda^- / [U{r, t) - [/(r, -oo)] d^r,(53) 
Js Jv 



where the expectation values refer to the ground state 
of the system consisting of the field, the atoms, and 
the reservoir oscillators responsible for the damping and 
Langevin forces acting on the atoms, and also 






The left-hand side of ([5^ is the electromagnetic energy 
fiux through a surface S enclosing the volume V, meaning 
that U{r,t) — U{r,—oo) is the energy density variation 
associated to this flux (through S) between i' = — oo 
and t' = t. The latter is immediately seen to vanish for 
static systems in equilibrium, but this can be circum- 
vented by imagining that we bring our system adiabati- 
cally from some reference configuration [2^ . We can then 
use the approximation that the system is stationary at 
both t' = — oo and t' = t, but as those instants are now 
characterized by different configurations, the subtraction 
in (|54p does not vanish anymore. In addition, by choos- 
ing our reference configuration as one where the different 
parts of the system do not interact, we can just subtract 
C/(r, — oo) from (|54p (as an irrelevant constant) and write 



C/(r,i) 



1 

8^ 






(55) 



We remark that the energy density [/(r, i), being a quan- 
tity that is associated with the electromagnetic energy 
flux, is not necessarily the same as the energy density 
u{y) of the previous section; in particular, as we shall 
see, it contains the energy of the bath field. In terms of 
the Fourier components of the fields [cf. ©], the total 
equilibrium energy density is 



"('-') ^s 



duj 



dio' ■ 



-i{uj'+Lj)t 



J UJ' +UJ 

X a;'(E(r,w) •D(r,a;')) +a;(D(r,a;) •E(r,w')) 
+ uj'{Ii{r,Lj) ■ B{r,uj'))+Lj{B{r,uj) ■ H(r, w'))] • (56) 

For zero temperature we can write this in terms of strictly 
positive frequencies as follows: 



I f°° f°° gH'' 
Ue{r,t) = — duj du' -- 

SttJo Jo ^ 



i{uj —cj)t 



c.(D(r,c.)-Et(r,c.')) 



(58) 



and, for the non-magnetic media under consideration. 



U^^{r,t) 



— / da;(H(r,L.).Ht(r,w)). (59) 



U{r,t)^Ue{r,t) + Um{r,t), 



(57) 



We have used the fact that E(r, — w) = E^(r,a;), etc. 
and, for uj,uj' > 0, 

(E(r,c.) • D(r,L.')) = (Et(r,c.) ■ t)\r,u;')) = 0, (60) 

and likewise for the corresponding bilinear magnetic field 
products. We have also made use of the fact that, for zero 
temperature, only the positive-frequency field operators 
Et(r, oj), etc. produce a nonvanishing result when act- 
ing on the initial state of the system, i.e., they act as 
creation operators, whereas the positive-frequency oper- 
ators E(r, cj) are annihilation operators and yield zero 
when acting on the initial zero-temperature state. 



B. Electric and Magnetic Fields 

The electric field operator E(r, w) satisfies 

2 2 

[V X V X -^e(r,w)]E(r,L.) = ^K(r,c.). (61) 
c c 

The solution for E(r, lo) in the medium can be written in 
terms of the Green dyadic satisfying (P7|) : 

Mr,Lo) = £^ydVG„(r,r',c.)K,(r',c.), (62) 

while for D(r, w) we have 

D(r,cj) = e(r,w)E(r,cj)+K(r,cj). (63) 

Based on the harmonic-oscillator reservoir model, it is 
shown in I, for instance, that the noise polarization K 
has the thermal equilibrium properties 

{K,{r,co)K,{r',co')) = {K}{r,co)K}{r' ,co')) = 0, (64) 
and 
{k}{r,u;)k,{r\u;')) = 4fie/(w)J,, J(l. - w')'5'(r - r') 



{K,{r,u;)kUr',u;')) = 4heiiLj)S,jSiuj - u;')S^{r - r') 



1 



^hw/ksT _ I 



(66) 



The expression for H(r, w) in the medium follows from 
the Maxwell equation V x E = -(l/c)9H/9t: 



H(r,w) = -i-V X E(r,a;) 

dV[VxG(r,r',a;)]-K(r',a;). 



(67) 



C. Total Energy Density 



Consider first the energy density Ue- The expression 
(p5)) for D(r,a;) results in two contributions to Ue- The 
contribution from e(r, aj)E(r, oj) is 

-1 ^oo (■oo i(ijj'—u)t 

Ui^\v,t) == —Re / dw / du'—, 

Stt Jq Jq w' - u; 

w'e* (w') - cje(a;)j (E(r, w) • Et (r, w')) 
-ReTr / dww^ / duj' 



Stt^c 



2„2 



(68) 



where to simplify notation we have suppressed the r de- 
pendence of e. We have employed the identity [l^ [see 
also Eq. (Ull) below] 

(E(r,a;)-E(r,a;')) - -^Tr G/(r, r,a;)5(tj - a;'), (69) 

with G/ denoting the imaginary part of the Green dyadic. 
In order to evaluate (p5)) we write the permittivity in 
terms of its real and imaginary parts, e(u;) = efl(w) + 
ie/(w)). In the term containing efl,(w) we use 



f{u:') - !{uj) _ dfiuj) 



lim 



du) 



(70) 



with /(w) = wefl(a;), while in the term containing e/(w) 
we use 



lim ■ = 2it. 



(71) 



After straightforward manipulations we obtain 
f^i'^W = ^Try dc.c.2£[c.eH(c.)]G,(r,r,c.) 

t /'OO 

+t—-rTr dwuj^ ei{iu)G I {r,r,ou). {72) 

The contribution to t/g (r, t) from the noise polarization 
part of D is 



U^'>{r,t) = Re- 



OO /"OO 



K2)^. ,^ _ R„^ /""rf^ r" dio'^ 







?'(ci;'— Ci;)i p 







Re — / dcj doj 

Stt Jq J^ lo' - uj An 

, V ,'2 



'(E(r,c.).Kt(r,w'))-w(K(r,a;)-Et(r,c^')) 



■' ' ^""^ /rfVG,,(r,r',a;)(i^,-(r,a.)i^/(r',c.')) 



2 , dVG^(r,r',w')(i^.(r,w)A-j(r',c.')) 



(73) 



r 



Using (j66p with T = 0, and the same sort of manipula- medium, and then 
tions used to evaluate C/i (r), we obtain 



f/P(r,i) 



? 9 



dujijp'ei iio) — — [wTr G^ (r, r, w)] 

OLlI 



'^^W = ^;^y dwc.2|_[^ej,(c.)]TrG/(r,r,c.) 
_e,(^)A[^TrGi,(r,r,L.)]}. (75) 



fc (.OO Integrating by parts the integral involving the second 

t / (icjo;' e/(w)TrG7(r,r, w).(74) term on the right, we obtain 

47r^c^ Jo 

h r°° 8 

+ 8^^ J dioujhi{uj)TrGR{r,r,u;). (76) 



Thus the terms proportional to t in ([7^ and (fM)) cancel, 
as discussed in I for the special case of a homogeneous 



10 



It follows similarly from ([59)) that 

h 



t/m(r) 



87r2c2 



Im / da;w2e(a;)TrG'(r,r,w). (77) 



The total zero-point energy density of the dispersive and 
dissipative dielectric medium is therefore 

X TrG(r,r,a;) 

h 



87r2c2 

n 

~ 87r2c2 JO 

+ ujv(i') 
u(r) + U7v(r), 



/•OO 

/ (ia;aj^e/(aj)TrGi^(r, r, oj) 
Jo 

|TrG(r,r) 



Im / dujLj' < — — lue] + e 



I du! 



(78) 



where u(r) is the energy density ([57]) obtained in Section 
II. 

Two "additional" contributions to the energy density 
have appeared naturally in the calculation of U{r) based 
on Poynting's theorem. The first one is the energy 
— i?absi [Eq. ([M)) ] associated with the absorptive loss 
of field energy, which is exactly cancelled by the energy 
+i?absi picked up by the reservoir oscillators, as already 
indicated in Eqs. d^jl-GSl) in Section [IlB. 

The second additional contribution, given by MAr(r), is 
exactly the same as in Eq. (PH]): 



WAr(r) 



/>oo 

■ / dwa;^e7(r,a;)TrGfl(r,r,w) 
Jo 



i — (E-K + K-E). 
2 87r ^ ' 



(79) 



As in the discussion that led to Eq. (|40p , this shows that 
WAr(r) is the polarization energy acquired by the "noise" 
dipole moments of the atoms as they interact with the 
electromagnetic field. 



D. Homogeneous Medium 

When spatial dispersion is negligible, equation ([75)) 
is a general expression for the zero-point energy den- 
sity of a linear dielectric medium, including the effects 
of the fluctuating electric and magnetic fields, the mate- 
rial medium, and the reservoir. We now show that this 
expression reduces to the expected energy density in the 
limit of a homogeneous medium. 

The (retarded) Green dyadic for a homogeneous dielec- 
tric medium is given by [3l| 



Gy(r,r',c^)= 1 



1 \ e 



ikfn |r — r'l 



(80) 



where A:^ — e(cj)a;^/c^. Thus, 
G,,(r,r',a.) = - 5.,,5^{v-v') 

Oft;™ 



RjRj/R — [dij — 3RiRj / R") 



x(l 



1 



where R= |r - r'|. For i? — ;• 0, 

47r 



aikmR 



R 



(81) 



G.j(r,r',( 



3fc: 



5..,5^{v~v') 



2 "^y 



1 



2i 

+ 3 "' '^ kim 



((5y - 3R^Rj/R^). (82) 



Ignoring for the moment the terms that diverge for i?, — >■ 
(see below), we put 



Giij{r,r,uj) = [2nR{u)uj/3c]6ij, 
Gffij(r,r,a;) = -[2n/(w)a;/3c](5y 

in ([75)) and obtain 

h 



(83) 



Z(r) 



dcjuj' 



■■{[2 



~ni — {uj€i) ~ 2njei\ 



^R + ^^R\nR 



(84) 



Evaluating the second term in curly brackets by partial 
integration we obtain 



u(r) 



27r2c 



dujLO ri^(a;)— — [ajri/j(a;)]. (85) 



This is exactly the energy density derived in I for the 
special case of a homogeneous dielectric medium. 

A few other points regarding the model of a homoge- 
neous medium may be worth noting. First, the Green 
dyadic ([8T)) for a homogeneous dielectric leads via equa- 
tion (|62p to the (Huttner-Barnett) expression (80) of I for 
the quantized (transverse) electric (and magnetic) field in 
a homogeneous dielectric, as is easily verified; these ex- 
pressions reduce to the standard ones in the limit of free 
space, as is shown by considering their limiting forms as 
tR — > 1,6/ — > 0. Second, equation ((53)) and the zero- 
temperature expectation value [cf. Eq. (82) of I] 



(E(r,c.).Et(r,a.')) = ^ 



nR{uj)5{uj ~ uj') (86) 



imply 



(E(r,a;)-Et(r,^')) = -^^G/(r,r,a;)5(c^-a;'), (87) 



TT C^ 



consistent with ([5^ . 

The delta function term in ([81]) and ([82)) has its origin 
in the singular divergence of R/R^ a.s R ^ 0, 



1 



w^-'.-. 



An 



3k 



S^,6^r-r'). 
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Physically, this is a direct consequence of the continuum 
and (spatial) local approximation, and can be remedied 
by allowing for spatial dispersion or/and the granular 
structure of matter [32[- As remarked earlier, such a sin- 
gularity is to be subtracted from G^. Its contribution 
to equation ([I} is a force density whose integral over all 
the space vanishes for bounded dielectric in a vacuum, 
as expected. Also, the \/B? divergence in ([5^ can be 
remedied by replacing the point dipole to which it cor- 
responds with a dielectric sphere of radius TZ and then 
allowing 7^—5-0; the result is |33t 



Gi.j{v,v',u) 



oh.2 ^J ^ 



r') 



(89) 



for |r — r'l infinitesimally small. Again, as noted ear- 
lier, this term is to be subtracted from Gij, and gives a 
vanishing contribution to the force density ([T]). 

Finally we note that throughout this paper we are 
working, albeit formally, with the full (transverse plus 
longitudinal) electric field and Green dyadic, as is re- 
quired when the permittivity can vary with r. If we con- 
sider, for example, the spontaneous emission rate A of 
a guest atom in a homogeneous dielectric, which is pro- 
portional to the imaginary part of G, we find that, in 
addition to the familiar contribution equal to the free- 
space emission rate times the real part of the refractive 
index at the transition frequency luq (and possibly also in- 
cluding local field corrections) , A has a contribution pro- 
portional to Im(l/fc^„) = (c^/ajQ)e7(ajo)/|e(wo)|^, where 
Wo is the transition frequency [cf. Eq. ((88)) ]. Follow- 
ing the approach of Barnett et al. [33[, this part of A, 
which is attributable to the longitudinal part of the elec- 
tric field, is found to depend on an unspecified parame- 
ter R related to the effective distance between the guest 
atom and the host medium. Without modifying the con- 
tinuum model to allow for such a parameter, we would 
obtain no such contribution to A. The same is true in 
I, where we considered only the transverse component of 
the field, in line with the assumption of a homogeneous, 
continuous medium made in that paper, where there was 
no allowance for the finite distances between atoms and 
the near-field interactions between atoms associated with 
longitudinal fields. 



E. Parallel Plates 

Equation ([T]) as such is not in general very useful for the 
calculation of Casimir forces between arbitrarily shaped 
media [SJ]. One exception is the force between two di- 
electric media occupying the half-spaces z < and z > d 
and separated by a third dielectric [ifl, UM ■ I^ ^^^^ ^^^^ 

e(r, uj) = ei{uj)0{~z) + e2(w)0(z - d) 

+ e3{uj)[e{z)~e{z-d)] (90) 

where 9{z) is the unit step function. Then the integral 
over space of the force density is simplified by the fact 



that Ve(r,a;) involves simple delta functions in z, yield- 
ing 

m) ^ I ^'r f (.) . -.^I,. [ ,u.^ coth [^) 

x{e3(w)G,,(0+,0+,w)-ei(w)G,,(0-,0-,w) 
+ t2{Lo)Gu{d+ ,d+ ,u:) ~ e^{w)Gu{d- ,d- ,Lo)Y (91) 

where Gi,;(0^, 0^,aj) and Gu{d^,d^,uj) stand for 
Gii{z -^ Q^,z' -)• 0=*=,^) and Gii{z -^ d^,z' -> d^,uj), 
respectively. The Green dyadic for this geometry can 
be evaluated without difficulty [1^ |35[ , and its trace at 
coincidence can be written as 



Gii{z,z,u)) 



I c 



u 



_ rd^ 

2tt oj^ei J Ki 

2 "1 



^l r n^i 



R\" 



'-^eiRf U-2»'^i^ (z < 0), 



(92) 
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27r w^es J K3 

o^tm tm 2iK3d tm 2iK3Z tm 2iK3(d-z) 

12 1 '2 



_2_ 



2 j,tc^tc 2iK,3d _|_ j,to 2iK3Z , ^tc2iK3 (d—z) 
2 ^tm^tm 2iK3d , ^tm 2iK3Z , ^tm ^2iK3(d-z) 

(0 < z < d), (93) 



Guiz,z,uj) 



K-^ ~r A^ii 



27r uj'^e2 J K2 

Je2i?2' [e^--^ (z > d), 



RT 



(94) 



where kj = . /ej-w^/c^ — /cy (the w-dependence of ej is 

implicit throughout), the Fresnel transverse electric (te) 
and transverse magnetic (tm) reflection coefficients are 
given by 

K3 — Kj ^j^3 ~ l^j 



and we defined 



-rp^r|e2-3d 
1 1 - rPrPe2'«3d ' 



m 



-rl + rle^^'^^'^ 
1 - r^rf e^^'^sd ' 



(96) 



(97) 
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Dp = 1 - rfr'^e^'^'"^ (p = tc, tm). (98) 

Substituting dMI-dMl) into dH]) we get, after straightfor- 
ward (if a bit long) manipulations 



/hA 1"°° 

cPr f (r) = -z ^"3 Im / duj cotli 

X / d^kii Hi ^ — 



2kBT 



— D„ 

p— tc.tni ^ 



(99) 



that is just the force between two parallel half-spaces 
separated by a dielectric [10| . 



IV. CASIMIR FORCE DENSITY FROM THE 
MAXWELL STRESS TENSOR 

The force on the material resulting from the Lorentz 
forces acting on the particles constituting the medium 
can be obtained from the Maxwell stress tensor [33] . We 
now show that this force is equal to the Casimir force 
density ([T]). 

Recall first the classical theory for the force density in 
a dielectric medium in which there are no free charges or 
currents. From the macroscopic Maxwell equations one 
derives a force density fi = djTij, where the stress tensor 
is 

^y = ^ l^^Dj + H,Hj - i (E • D + H . H),5,,] . (100) 



It follows that 



h^—[{d,E^)D^-Ej{d,Dj)] 
an 



(101) 



In identifying f (r) as the force per unit volume on the 
medium at the point r we are assuming conditions such 
that (9/9t)(D X H)/47rc can be taken to be zero, since 
we are dealing with a stationary equilibrium situation 
(see also [3^). We continue to restrict ourselves here, 
as in I, to isotropic media, in which case the Minkowski 
form of the stress tensor defined by pOO|) is symmetric, 
consistent with angular momentum conservation. 

In quantum theory we replace Ej and Dj in poip by 
operators, symmetrize, and take expectation values: 



fi == _Rc 

•^ Sir 



{{d,E,)D,)~{E,{d,D,)) 



(102) 



or 



/.(r) = —Re 

OTT 



duj / duJ'U[^^Ej{v,L0)]bJ{v,L0')) 

{E,{r,u;mD,{r,u;')])}. (103) 



Since the electric displacement vector D(r,w) = 
e(r,w)E(r,w)-|-K(r,w), 

{E,{r,cj)D,{r',Lo')) = e{r' ,Lo'){E,{r,u;)E,{r' ,u;')) 

+ {E,{r,uj)K,ir',cv')). (104) 



In the second term on the right wc use (|62|) . while in the 
first term we use (1691). Then 



Rc{E,{r,uj)I),ir',u;')) = — Im[6(r',^)G,,(r, r',c^)] 

x5{lu-lo'), (105) 

which leads us again to the force density ([T|) : 
/*(!■) = -^Ini / duj duj' 

J^TT Jo Jo 

X lim {{[d,E,{r,Lo)]D,{r',u')) 
-{E,{r',u)[d,D,{r,u')])} 

/•CJO 

^Im / diouj'^die{r,Lu)Gjj{r,r,Lo). 
Jo 



(106) 



V. CONCLUDING REMARKS 

In Section |ll] we considered a collection of station- 
ary electrically polarizable particles in free space and 
at zero temperature and obtained formal expressions for 
the many-body van der Waals interaction energies among 
these particles. The results are equivalent to those ob- 
tained by various other methods, but unlike most of ear- 
lier work we accounted explicitly for dissipation by cou- 
pling each particle to a reservoir of harmonic oscillators 
rather than by invoking analytic properties of the permit- 
tivity. Within the continuum approximation we obtained 
expressions for the electromagnetic energy density and a 
Casimir force density; as shown in the Appendix, these 
expressions are easily generalized to the case of thermal 
equilibrium. The force density obtained is in agreement 
with that originally derived by Dzyaloshinskii et al. jlOf . 
In Section Hill wc derived the formula ([57)1 for the energy 
density starting from the continuum approximation and 
Poynting's theorem in macroscopic quantum electrody- 
namics, and in Section HVl wc derived the formula ([1]) for 
the force density using the Minkowski form of the stress 
tensor. 

Under the assumption of a continuous medium, and ig- 
norin g sin gle-particle self-energies, the expression (|571) is 
exact |38| . Similar expressions are well known in classical 
electromagnetic theory 39] as approximations applicable 
for frequencies at which absorption is negligible. The rea- 
son our expression is exact is not, strictly speaking, due 
to our quantum treatment of the field, but rather to the 
fact that Fourier components of the zero-point field field 
at different frequencies are uncorrelated |8| (as they are 
at finite equilibrium temperatures). 

Under the (unrealistic) assumption that absorption 
can be completely ignored, so that e is real, we would 
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write 
u{r) 



1 

8^ 



r d -I 



x(E(r,w)-Et(r,,^)) 



-Im 



d 



87r2c2 

X Gji(r,r,a; 



I dwa;2| — [we(r,a;)]+e(r,w)| 

(107) 



Since both e and G as functions of complex frequency 
are analytic in the upper half of the complex frequency 
plane, we can perform a Wick rotation and write 



u(r) 



A^lmJ^^ d^e{§-^Mr,^0]+e{r,^0} 



8 



X Gu{r,r,iO, 



(108) 



where wc have used the relation between (E(r,a;) • 
E^(r, w)) and G{r,r,Lu) for purely imaginary frequencies 
i^. In other words, consistent with the remarks in the 
Introduction, we can obtain the correct energy density 
from the formula (|107p . which assumes there is no ab- 
sorption, by formally replacing e(r, co) and G(r, r, w) by 
their (real) values e(r,i^) and G(r, r, i^). Similarly we 
can formally interpret the energy density p7p as an in- 
tegral over zero-point energies {l/2)hh' by writing it as 

u{r)=J di^p{r,iy)(^^hi^y (109) 



where 



is a local density of states [41 



'^{A.[^,^r,^,.)]+eir,^,.)]} (110) 



proach of Scction|TT]wc replace (E(r) • E'l'(r)) in Eq. ([14)) 
by (E(r) • Et(r) + Et(r) • E(r)), since the quantized field 
E involves photon annihilation operators and docs not 
give zero when acting on a thermal state. We similarly 
replace (fTS]) bv figf 



(Eo(r,L.) • Ej(r,w') + El{r,Lu) ■ Eo(r,c^')) = 



— coth 



2kBT 



')lmTrGo(r,r,a;)(5(w-a;'). (HI) 



Then formulas such as p7)) - ((T9| and ([37]) are general- 
ized to finite temperature by simply inserting the factor 
cot]i{huj /2kBT) in the integrations over lo. Thus, for ex- 
ample, the energy density (P7)) generalizes to 



M(r) 



87r2c2 



ImTr / dujuj'^ coth ( — ) 

\2kBTJ 



X [2e(r,a;) + a;e'(r,w)]G(r,r,a;). (112) 



The same arguments leading from the energy density to 
the force density apply regardless of the temperature, so 
that the generalization of ([1]) is simply 

f(r) = -^Im^°°d..^coth(^)ve(r,.) 

xG,,(r,r,w). (113) 

We can convert in the standard fashion the integral over 
the real lo axis to an integral over the imaginary uj axis. 
Using the fact that e(r, w) and G(r, r, w) are analytic in 
the upper half of the complex frequency plane, and that 
coih{huj/2kBT) has poles at ci;„ — {27TikBT/h)n = i^„ 
for all integers n, the result is that the integration over 
w is replaced by a sum over n: 
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Appendix. Generalization to Finite Temperature 

It is an easy matter to generalize our approach and 
results to finite equilibrium temperatures T. In the ap- 



dujf{uj) 



2'KkBT 



n=0 



li^u), 



(114) 



where /(w) denotes the integrand in (|113p and the prime 
on the summation sign indicates that the n = term is 
multiplied by 1/2. Then 



f(r) 



kBT 
47rc2 



^a;^G,i(r,r,ia;„)Ve(r,ia;„), (115) 



n=0 



which is the result of Dzyaloshinskii et al. quoted earlier 
[Eq. P^]) ] . The same finite-temperature results are easily 
shown to emerge from the approach of Section IIIII when 
the temperature-dependent terms in the noise polariza- 
tion correlation functions (pS)) and (|M|) are retained. 
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